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Abstract. The main goal of the Andyrchy-BUST 
experiment is to study the primary cosmic rays 
spectrum and composition around the knee. The 
experimental data on the knee, as observed in the 
electromagnetic and high energy muon components, 
are presented. The electromagnetic component in our 
experiment is measured using the "Andyrchy" EAS 
array. High energy muon component (with 230 GeV 
threshold energy of muons) is measured using the 
Baksan Underground Scintillation Telescope (BUST). 
The location of the "Andyrchy" right above the 
BUST gives us a possibility for simultaneous mea- 
surements of both EAS components. 

Keywords: EAS, primary spectrum, primary com- 
position 

I. Introduction 

In the range of primary energies of 10 14 — 10 15 eV per 
nucleus, direct methods for studying the energy spectrum 
and nuclear composition of primary cosmic rays become 
inefficient because of a decrease in the flux of primary 
particles with an increase in their energy. Therefore, at 
these and, of course, higher energies, indirect methods 
based on simultaneous measurement of the characteris- 
tics of different components of extensive air showers 
(EASs), which are initiated by the primary particle 
in the atmosphere, are used. But the interpretation of 
these measurements requires their comparison with EAS 
simulations in the atmosphere. In turn, the calculation 
results depend on the hadronic interaction models. The 
main problem is the extrapolation of these models 
into kinematical and energy regions still unexplored by 
present-day collider experiments. So, the measurements 
of different EAS components are now used for both 
studying the primary composition and testing interaction 
models Q - Q. 

In this paper three types of experimental data 
are analyzed: muon number spectrum, EAS size 
spectrum and correlation between muon number and 
EAS size simultaneously measured. Integral muon 
number spectrum has been measured using the Baksan 
Underground Scintillation Telescope (BUST) (8). The 
EAS size spectrum has been measured using the 
"Andyrchy" EAS array (9J. The dependence of the 
mean number of high energy muons on EAS size 
has been measured by simultaneous operation of both 
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Fig. 1: Andyrchy EAS array. 



devices ifTOl. ifTTI. 



II. Facilities 



The "Andyrchy" EAS array is located on the slope 
of the Andyrchy mountain above BUST (43.3° N, 42.7° 
E) and consists of 37 plastic scintillation detectors. A 
plastic scintillator of the 5 cm thickness has an area of 
lm 2 and is viewed with a single PMT. The detectors are 
designed for both timing measurements (for EAS arrival 
direction) and evaluation of primary energy (via EAS 
core localization and determination of total number of 
particles in shower). The distance between the detectors 
is about 40 m in projection to the horizontal plane and 
the overall area of the installation is 5 • 10 4 m 2 ; the 
solid angle at which the array is viewed from the BUST 
is 0.35 sr. The central detectorof the array is above the 
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• simulated spectrum: FfN^) ~ N 
X reconstructed spectrum 




Fig. 2: Simulated and reconstructed size spectra. 



BUST'S center at a vertical distance of about 360 m and 
at 2060 m above the sea level. The difference between 
the heights of the upper and lower rows of detectors is 
-150 m. 

The energy deposition measurement is performed in nat- 
ural units, so called relativistic particles. One relativistic 
particle (r.p.) is the most probable energy deposition 
from a single cosmic ray particle. For our detector it 
is 10.6 Mev |[T2l . The range of the energy deposition 
measurements is from 0.5 r.p. (the threshold of the 
Charge-to-Time Converter) up to more than 1000 r.p. 
Trigger formation and all measurements are performed 
in a registration room, which is placed near the center 
of the array (length of connection cables is up to 280 
m). The shower trigger condition requires signals from 
4 detectors within 3 microseconds. The trigger's rate 
is about 9 s _1 . The array and its characteristics are 
described in more details in iTTZll . 

BUST lfl3l is a large device 16.7 x 16.7 m 2 area 
and 11.1 m height), located in a cave under mountain 
slope. The four vertical sides and four horizontal planes 
are completely covered with standard liquid scintillation 
detectors. The standard detector consists of an alu- 
minium tank with 0.7 x 0.7 x 0.3 m 3 dimensions and is 
filled with liquid scintillator on the base of white-spirit. 
Total number of the detectors is 3180. Every counter 
is viewed with one PMT. The construction of BUST 
allows one to reconstruct tracks of muons crossing the 
telescope. Coordinates of hit detectors is used as input 
information for muon group parameters determination. 
The telescope allows one to determine the number of 
passing muons, their coordinates (with 0.7 m accuracy) 
and the arrival direction (with 1.5 degree accuracy). 
The coincidence trigger between "Andyrchy" and BUST 
is produced when one or more muons crossing the 
telescope (« 12s -1 ) coincide with the shower trigger 
within 51.2 microseconds; the coincidence rate is about 
0.1s- 1 . 



Fig. 3: EAS size spectrum. Points - experiment. Lines 1 
and 2 are calculated spectra for the primary compositions 
1 and 2, correspondingly. 

III. EAS SIZE SPECTRUM 

The standard definition of the shower size N e is 
the total number of the charged particles (mainly e^) 
at the level of observations. As a scintillation detector 
measures the energy deposition, and not the number 
of particles, the reconstruction of shower parameters is 
performed in units of relativistic particles (r.p.). The 
measured size iV r p . is the total energy deposition in 
allegedly continuous infinite detector. The shower size 
N r , p ,, the slope of the lateral distribution function and 
the core location are determined by a x 2 -hke method, 
in which the logarithm of the energy deposition in each 
detector is compared with the one expected from the 
NKG lateral distribution function 

with ro = 96 m. The NKG function reproduces with a 
good accuracy the experimental data (9j. 

In the present analysis, only the showers with: 

1) sec 9 < 1.05 (near vertical events); 

2) > 22 fired detectors; 

3) > 4 detectors with energy deposition Ed > 10 r.p. 
well inside the array; 

4) reconstructed axes in central part of the array (the 
distance from the center is not larger than 50 m) 
were taken into account. 

The accuracy of reconstruction was calculated using 
data obtained from a simulation that includes the ex- 
perimental dispersion. Figure [2] shows the reconstructed 
size spectrum together with the simulated one, the size 
spectrum is reconstructed without distortions for showers 
with 5.75 < IgAV.p. < 8.0. For these showers the 
accuracy of the N r , p , determination is better than 15% 
and the accuracy of the axis position determination is 
better than 5 m. Figure [3] shows the measured differential 
size spectrum taken during live time 1.297-10 8 s (1501.2 
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Fig. 4: Integral muon tracks number spectrum. Points - 
experiment. Lines 1 and 2 are calculated spectra for the 
primary compositions 1 and 2, correspondingly. 

days), in r.p. units. The steepening of the spectrum is 
observed at lg N r . p . « 6.35. 

IV. Muon number spectrum 

Coordinates of hit BUST detectors are used to recon- 
struct tracks of muons crossing the telescope. Generally 
the number of muon tracks m differs from the number 
of muons in the group. In the case where the 
distance between muons is small enough (compared 
to the individual detector size) the number of recon- 
structed muon tracks is smaller than the number of 
muons in the group. Opposite effect is also present 
since interacting muons can produce particles which 
might increase the number of hit detectors, therefore 
the number of reconstructed muon tracks in such a 
case can be larger than the number of muons in the 
group. Furthermore there is some arbitrariness in muon 
track determination: for example, a track may cross two, 
three or four telescope planes and so on. Hence it is 
necessary to convert the number of reconstructed muon 
tracks to the number of muons in the group taking 
into account all mentioned effects. Conversion factors 
depend on muon lateral distribution function. Lateral 
distribution might depends on muon energy distribution 
in EAS and on primary nucleus energy (per nucleon) 
etc. In order to avoid additional uncertainty we use only 
experimental muon tracks number spectrum for study of 
primary composition. The conversion of the number of 
muons to the number of reconstructed muon tracks is 
included in the calculations. The integral spectrum on 
the number of muon tracks for near vertical directions 
(8 < 20°, effective muon threshold energy is 230 GeV) 
was measured for m = 1 - 250 ((Fig. |4). This spectrum 
was obtained using two BUST data sets. The first data 
set has been taken during 2001 - 2004 (live time 3.3 
years) and contains information about all the BUST 
events. The second data set has been taken during 1984 
- 1995 (live time 9.8 years) and contains information 
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Fig. 5: Nfj,(N rp ) dependence. Points - experiment. Solid 
line 1 and dashed line 2 are calculated dependencies for 
the primary compositions 1 and 2, correspondingly. The 
dotted lines are calculated dependencies for pure protons 
and iron nuclei. 



about BUST events where 100 and more detectors were 
fired. The later condition corresponds to 20 or more 
muons crossing BUST for near vertical directions. It is 
to be noted that owing to muon number tracks saturation 
there is no sense to analyze this spectrum for m > 170, 
where the relation m^/m becomes more than 1.5. 

V. The mean number of muons vs. EAS size 

The size, axis position and the EAS arrival direction 
are determined using the "Andyrchy" array data; the 
BUST data are used to determine the number of muons 
crossing BUST. The underground telescope measures 
only a part of the total number of muons in EAS and 
the uncertainty in the determination of the EAS axis 
position at the BUST level is comparable with the size 
of the BUST. The mean number of muons in BUST is 
determined as follows. Events within a given range of 
N e are grouped according to the distance to the BUST'S 
center with step SR = 10 m. For each group, the number 
of muons is 

Ki 

where K{ — K(Ri) is the number of EAS in the i- 
th group, and m,j is the total number of muons in the 
BUST in the z-th group for the j-th EAS. 

Thus, the mean number of muons for a given range 
of N r _ p _ can be written as 

M[Rj) 1 ^ 

1 1 3=1 

The mean number of muons in a shower is then calcu- 
lated as 
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TABLE I: Primary composition 1. 



z 


1 


2 


6-8 


10-16 


20-26 


~z 


1 


2 


7.2 


12.7 


25.2 


A 


1 


4 


14.4 


25.5 


54.2 


F z °, (m 2 s sr TeV)' 1 


0.12 


0.056 


0.03 


0.035 


0.0267 


1z 


2.75 


2.64 


2.66 


2.70 


2.62 



N, 



(3) 



where St = 200 m 2 is the effective area of the telescope, 
S r (Ri) is the area of the ring of radius R4 with SR = 10 
m. 

The fraction of muons in the telescope A(R) was 
measured as a function of the distance R between the 
center of the telescope and the EAS axis for a set of 
showers with a given value of N r . p , . This fraction for a 
given range of N r , p . is defined by 



n(R) 



(4) 



where n(R) is the mean number of muons in the 
telescope at the distance R, N^Nr.p.) is the mean 
number of muons in EAS. Because A(i?) depends on the 
lateral distribution function (LDF) of high energy muons 
and BUST geometry only, LDF can be recovered from 
these measurements. Preliminary results on high energy 
muons' LDF can be found in lfl4l. 



VI. Calculations 

The development of EAS in the Earth's atmosphere 
have been simulated by means of the CORSIKA code 
(version 6900) fl3]. The QGSJetII-03 and Fluka were 
used as the high and low energy hadronic interaction 
models. The CORSIKA output files were used then as 
input files for AndyrDet code, which performs a detector 
response simulation. The results of the simulations were 
summarized as parametrization functions of the EAS 
characteristics, that then was used for calculations of the 
observables (integral muon number spectrum, EAS size 
spectrum, N fl (N rp ) dependence). Vertically arriving CR 
particles were used as primary CR particles. Simula- 
tions were performed for nuclei with atomic number 
A = 1, 4, 14, 28, 56. Primary energy per particle Eq was 
taken from the range 10 4 GeV - 10 7 5 GeV with a step 
of 0.5 of energy decade. 

The integral muon number spectrum in BUST was 
calculated numerically in the same way as in (8). This 
calculation method needs only such characteristics of 
high energy muon component of EAS as: 1) N^(Eq,A) 
- muon production function (MPF) or mean number 
of muons per EAS produced by nucleus with atomic 
number A and primary energy Eq; 2) f(r,Eo,A) - 
lateral distribution function (LDF); 3) G(A,N fJl ,N fJl ) - 
fluctuation function (FF). 



p 

Fe 



log10(N 



Fig. 6: The number of high energy muons (with 
E m u > 230GeU) as function of N rp for the different 
primary energies of the primary protons and iron nuclei. 



The EAS size spectrum can be presented as: 

00 

dF A (Eo). 



i(N r . p .) = Yl [ ■ 

A J 



dEn 



-W A (Eo, N r . p .)dE (5) 



where dFA(Eo)/dEo is the energy spectrum of the 
primary nuclei and Wa(Eo, N r . p .) is the probability for 
primary nucleus with energy Eq to produce EAS with 
size N r , p , at observation level. 

The dependence of the mean number of high energy 
muons on EAS size was calculated taken into account 
the energy spectra of the primaries and anticorrelation 
between the number of high energy muons in EAS and 
EAS size at fixed primary energy (Fig. [6]). 

The calculations of the observables were performed 
for two composition models with five groups of primary 
nuclei. Energy spectrum of every primary group is a 
power law with rigidity dependent knee E^ z = Ek p ■ Z: 



dF z (E ) 
dE 



F z ■ Eq 



1 



Ekz 



(6) 



where Eq is energy per particle and F z ° is absolute 
flux at 1 TeV per particle ifTBI . ifTTl . The sharp knee 
was applied for both composition models: e c — 3.5 and 
7 C = 5.2. Both composition models do not give closest 
fit for the complete data set of direct measurements, 
but the data do not contradict these models. All-particle 
energy spectra for these composition are shown in Fig. 
[7] together with averaged spectrum from air shower 
experiments ifTTll . 



PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



5 



TABLE II: Primary composition 2. 



z 


1 


2 


6-8 


10-16 


20-26 


~z 


1 


2 


7.2 


12.7 


25.2 


A 


1 


4 


14.4 


25.5 


54.2 


F z °, (m 2 s sr TeV)' 1 


0.102 


0.079 


0.018 


0.018 


0.02 


72 


2.77 


2.65 


2.70 


2.70 


2.64 



10' t, 




log10(E /A) 



Fig. 7: All-particle energy spectra: lines 1 and 2 - for 
compositions described in Tables 1 and 2, correspond- 
ingly. Points - averaged all-particle energy spectrum 
from air shower experiments ifTTl . 



Fig. 8: Muon production function. Points - CORSIKA 
results. Line 1 - fit of the CORSIKA results. Line 2 - 
changed MPF 



The first composition is presented in Table 
1, for this composition total flux at 1 TeV is 
0.2677 (to 2 s sr TeV)^ 1 and protons knee position 
is Ek p — 4 • 10 3 TeV. This composition gives a good 
fit for the EAS size spectrum (Fig. [3] line 1). But the 
muon number spectrum calculated for this composition 
is in contradiction with experiment (Fig. |U line 1). 
For the second composition (Table 2) the protons knee 
position is E kp = 2 • 10 3 TeV and total flux at 1 TeV 
is 0.215 (to 2 s sr TeV)^ 1 . This composition gives a 
satisfactory fit for the muon number spectrum (Fig. |U 
line 2). But the EAS size spectrum calculated for this 
composition is in contradiction with experiment (Fig. [3] 
line 2). For both compositions the calculated iV M (jV r p) 
dependences are very close to one another and both are 
in contradiction with experiment (Fig. [5]). 

So, predictions from both considered mass composi- 
iton models do not give satisfactory fits to all data set 
obtained from our measurements (see Fig. [3]and|4]i. 

VII. High energy muon production function 

Discrepancies between EAS size spectrum and muon 
data can be lessened to a considerable degree if we use 
another muon production function for muons with E p > 
230 GeV . The MPF can be expressed by next formula: 

Fit of the CORSIKA (QGSJetII-03) results was obtained 
for parameters: b = 0.0018, c = 13.5, a = 0.43 and 
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Fig. 9: Integral muon tracks number spectrum. Points 
- experiment. Lines - calculated spectrum for the first 
primary compositions with: 1 - MPF from CORSIKA 
QGSJetII-03; 2 - changed MPF. 



j3 = 1.675 (line 1 in Fig. [8]). It gives asymptotic behavior 
for N„ ~ £ - 720 . 

The agreement between EAS size spectrum and muon 
data can be obtained for MPF with parameters: b = 
0.0035, c = 14, a = 0.42 and /? = 1.54 (line 2 in 
Fig. [HJ. The asymptotic behavior of the mean number 
of muons for this MPF is: ~ Eg- 647 . The integral 
muon number spectrum calculated for the first primary 
composition and both MPF's is shown in Fig. [9] 

The N ll (N rp ) dependences calculated using changed 
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MPF are in Fig. [10] One can see that the dependence 
calculated for the first primary composition and changed 
MPF gives acceptable fit with experimental data. 

VIII. Conclusion 

It is widely known that none of the present interaction 
models can completely describe a full set of experimen- 
tal data for cosmic rays. Joint analysis of the various 
characteristics of different EAS components, especially 
measured in one and the same experiment, can be used 
for both studying the primary composition and testing 
interaction models. 

In this paper three types of experimental data, taken in 
our experiment, were analyzed: high energy (E^ > 230 
GeV) muon number spectrum, EAS size spectrum and 
dependence of the mean number of high energy muons 
on EAS size. CORSIKA code v.6900, with QGSJetll- 
03 and Fluka as the high and low energy hadronic 
interaction models, has been used for EAS simulations. 
As it was mentioned above, the calculations of these 
observables need only the following characteristics of 
EAS: 

1) Wa(Eq, N r . p ., Nfj) - the probability for primary 
nucleus with atomic number A and energy Eq to produce 
EAS with size N r . p . and total number of high energy 
muons at the observation level; 

2) Wa(E , N r . p .) - the probability to produce EAS with 
size N r , p . at the observation level; 

3) N^Eo, A) - high energy MPF; 

4) G(A, NftjN^) - high energy muons fluctuation func- 
tion; 

5) f(r, Eq, A) - high energy muons LDF 

It should be noted that characteristics 2-4 can be 

derived from 1. 

The analysis has shown that: 

1) our experimental data can be brought into good 
enough agreement using changed MPF (Fig. [HJ; 

2) in both cases the primary composition gets heavier 
across the knee. 
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